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Abstract 

The Jahn-Teher distortive transition of LaMnOs is described by a modihed 3-state Potts model. 
The interactions between the three possible orbits depends both on the orbits and their relative 
orientation on the lattice. Values of the two exchange parameters which are chosen to give the 
correct low temperature phase and the correct value for the transition temperature are shown to 
be consistent with microscopy theory. The model predicts a first order transitions and also a value 
for the entropy above the transition in good agreement with experiment. The theory with the same 
parameters also predicts the temperature dependence of the order parameter of orbital ordering 
agreeing well with published experimental results. Finally, the type of the transition is shown to 
be close to one of the most disordered phases of the generalised Potts model. The short range 
order found experimentally above the transition is investigated by this model. 
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I. INTRODUCTION 



Orbital ordering and Jahn- Teller ( JT) distortions are observed in some perovskite-type 
3d transition-metal compounds such as the vanadates (e.g. V2O3/ LiVOs,— and LaVOs.— ) 
cuprates^ (e.g. KCuFs) and derivatives of the colossal magnetoresistive in the compound 
LaMnOa ^ which is the subject of this paper. 

A single electron, or hole, in cubic symmetry will be in a two-fold or three-fold degen- 
erate state if it occupies a eg or t2g orbit. This degeneracy will be lifted at a transition where 
the local symmetry is lowered and long range orbital order occurs in the low temperature- 
phase^. The local orbital order will be accompanied by local structure distortions because 
of a linear coupling between the occupations of the different orbital states and the lattice^. 
This is the JT effect. In a perovskite structure this is a distortion and rotation of the oxygen 
octahedra where in the lattice structure ABO3 (where A= rare-earth cation and B=Mn) the 
B atom is octahedrally coordinated to the oxygen. The local distortion around one B ion 
affects other ions in the crystal because the oxygen ions are shared between two octahedra. 
If a large fraction of the energy of the orbital ordering comes from the energy lowering due to 
the interaction with the lattice the interaction is known as a cooperative JT effeclf^^. There 
is a linear coupling between an eg orbital and the E distortion modes of the octahedron. 
The coupling may be determined from a classical harmonic approximation to the lattice 
dynamicsifl'. 

The subject has attracted much interest recently because of the interplay between 
the magnetism and the lattice. The coulomb repulsion alone which gives rise to magnetic 
superexchange^ can also drive orbital order and hence there has been much efforts to de- 
termine the relative importance of the lattice effects, cooperative JT, and the electronic 
energies in driving the transition. 

We consider LaMnOa , where Mn^+ has a c?^ configuration, i.e., four electrons in d 
orbitals. The wave functions, d^^-yZ and d^^'i-r^ , are called eg orbitals, whereas, dxy, dyz, 
and dzx are called t2g orbitals. The level splitting is shown in Fig. ^ in a cubic crystal and 
with a tetragonal distortion. The orbital ordering is shown in figure ^p. When electrons 
occupy these wave functions, the ground state is determined by the semiempirical Hund's 
rule in LaMnOs because the exchange energy is greater than the crystal field energy, 6cf- 
Hund's rule coupling causes all of the spins of the electrons to be aligned parallel, that is 
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(a) 



FIG. 1: (a) Splitting of 3d-levels in a cubic crystal field with an additional tetragonal distortion. 
Electron occupation of the four 3d-electrons in Mn^+ is shown by arrows, (b) the orbital ordering 
of the eg states in an octahedral crystal of LaMnOs at T < Tjt- One of the two bonds of Mn-O-Mn 
bond in the x — y plane contracts and the other expands. The La'^^ ion in the centre of the cell is 
omitted for clarity 

5 = 2, where three electrons are in the t2g orbitals and one electron occupies one of the Cg 
orbitals. The MnOe octahedra in the low temperature phase have two short, two medium, 
and two long bonds^i while above Tjt the octahedra have six almost equal bond lengths. 

This paper develops a model for the temperature dependence of the orbital order- 
ing. There have been several studies of the effect of temperature in the orbitally ordered 
phase. The orbital ground state is well known as a C-type, antiferromagnetic in x-y plane 
and ferromagnetic in z-direction, with long-range ordered, JT distorted, MnOe octahedra^*^ 
as sketched in Fig. [T)d. The elongated occupied Cg orbitals alternate between pointing 
along X and y directions, the so-called orthorhombic structural phaso^i^. Above Tjt, the 
structural phase is pseudocubic with almost regular MnOg octahedra. Also above Tjt the 
long-range orbital ordering disappears. The structure was studied using the thermal analy- 
sis high-resolution neutron-powder diffraction patterns^ and Raman phonons as function of 
temperature^^. Evidence has also been presented^ for orbital ordering in LaMnOs which is 
based on the splitting of the Mn 4p levels by the Mn 3d orbital ordering. 

LaMnOs undergoes a structural phase transition^ii^ at the JT temperature, Tjt = 
750K. Murakami et al.^ showed that the order parameter of orbital ordering as a function 
of temperature decreases rapidly for T > QATjt and disappears at Tjt = 780K in a first 
order transition. These results have been re- interpreted and fittedi^ by a critical exponent. 
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(3 = 0.16(1), at TjT when the temperature dependence of the electron-spin resonant in 
Lao.gsSro.osMnOs has been investigated. 

Sanchez et al^ showed by means of X-ray absorption near edge structure and extended 
X-ray absorption fine structure at the Mn i^'-edge that the structural transition originates 
by the ordering of tetragonally distorted octahedra and they analysed the entropy content 
of the transition within the framework of the three-state Potts modeU^*i2ii§ with nearest- 
neighbour antiferrodistortive coupling. It is confirmed^S,, by using neutron diffraction data 
and a combined Rietveld and high real space resolution atomic pair distribution function 
analysis, that the nature of the JT transition around Tjt is to be orbital order-disorder and 
the intermediate structure suggests the presence of local ordered clusters implying strong 
nearest-neighbour JT antidistortive coupling. 

Millisi^ derived a classical model, which was based on previous work by Kanamori^S, 
for the lattice distortions in manganites. The model may be approximated either by an anti- 
ferromagnetic xy model with a modest threefold anisotropy or by a 3-state Potts model with 
an antiferromagnetic first-neighbour interaction and a weak second-neighbour interaction. 
This differs from our model which deals only with the nearest neighbour interaction. 

We study the orbital ordering appropriate for LaMnOs obtained from our model and 
the JT distortion which occurs above JT temperature in that phase. The model proposed 
here for the disordered phase is an array of localised distortions with no long range or- 
der. This differs from the model of Zhou-Goodenougb^i who postulated that dynamic JT 
effects occur in the high temperature phase. Zhou and Goodenough suggested that the 
orbital order- disorder transition at Tjt is to a dynamic JT stabilisation at the Cg electrons 
persisting into the high-temperature 0-orthorhombic phase, and they proposed a vibronic 
mechanism to explain the electrical resistivity and thermoelectric power above and below 
Tjt respectively. 

It is well known that the phenomenological models such as pseudo spin or Potts models 
are very useful to study order-disorder transitions. In the LaMnOs system, the MnOe octa- 
hedra are not independent as they share oxygen atoms with their nearest neighbours. The 
standard 3-state Potts model alone does not have the correct orbital ordering for LaMnOs 
as a possible ground state. In this paper, we describe a phenomenological model that gives 
the correct ground state for the orbital ordering and a good value for high-T entropy for the 
LaMnOa phase. The model is called the anisotropic Potts model and has been presented 
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earlier—. 

In section ini the anisotropic Potts model and its phase diagram are presented and its 
application to LaMnOa is in section 11111 The model is used to give understanding of the 
orbital ordering just above Tj^ in section IIVI Finally, we conclude the results in the last 
section. 

II. METHODOLOGY 
A. The model 

We model the transition at Tjt considering only the largest distortions that leave the 
crystal tetragonal in the low temperature phase. If the tetragonal axis is chosen along z 
direction the orbital ordering is characterised by the vector Q = (110). There is a staggered 
order in the x — y plane and the orbits are stacked ferromagnetically up the z axis as shown 
in Fig. ^b). However, in Fig. ^b) we showed the ordering for a pure JT system where the 
orbits shown were |a;) = |3x^ — r^) and \y) = |3?/^ — r^). These orbits are found from the Cg 
doublet for the case 9 = 

\9) = cos-\3z^ — r^) + sm-|x^ — y^). (1) 

Experimentally^ it is found that if the tetragonal axis is along z direction, then, the stag- 
gered orbits are given by | ± Oexp) where d^xp is closer to 7r/2 than 27r/3. We note that, 

I ± — ) = cos—\3z — r ) + szn—\x — y )■ (2) 

This gives, when | — |) is used, ^|2.7x^ — 0.7?/^ — 2z^) which has its main lobe along x 
direction. Similarly, using | + |) gives a state, ^|2.7?/^ — 0.7a;^ — 2z^), which has its main 
lobe along y direction. We use the JT states in our phenomenological model because these 
allow us to use the same states below Tjt for any of the three possible Q vectors, (110), 
(101) and (Oil), and also above Tjt- The model presents the correct dimension of the order 
parameter. 

There are six equivalent orderings that can occur. These correspond to the three 
choices for the tetragonal axis and, then, the phase (±1) of the order parameter. In the 
high temperature phase of the Potts model the orbits are in one of states 9 = {] oi 9 = 
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This differs from Zhou and Goodenougb^i who assumed that above Tjt the orbit is rotating, 
thus, all values of 6 are accessed dynamically. Because the orbital ordering in the transition 
metal oxides has three anisotropic states it is natural to set up a three states Potts model. 
The standard g-state Potts model-*^^ consists of a lattice of spins , which can take q different 
values from 1 to q, and whose Hamiltonian is 



J ^ 

(id) 



(3) 



where Si = 1,2, is one of the q states on site i, Ss^^Sj is the Kronecker function which 

is equal to 1 when the states on sites i and j are identical. Si = Sj, and is zero otherwise, 
means that the sum is over the nearest neighbour pairs, J is the exchange integral 
and is the total number of sites in the lattice. For q = 2, this is equivalent to the Ising 
model. The Potts model is, thus, a simple extension of the Ising model, however, it has a 
much richer phase structure, which makes it an important testing ground for new theories 
and algorithms in the study of critical phenomenai^. 
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(b) Ji (disfavoured) 



(c) Js (favoured) 



FIG. 2: The exchange interactions used in the anisotropic Potts model to give orbital ordering in 
LaMnOa , (a) if a Mni-0 bond expands, the oxygen shared 0-Mn2 bond should contract, (b) Two 
identical orbits in the same lattice vector, Ji > 0, disfavoured, (c) Two identical orbits in different 
lattice vector, J2 = 0, favoured. 

We have presented^^ a modified Potts model on a cubic lattice where the interaction 
between two sites depends on both the orbitals on these sites and also the directions of the 
nearest neighbour bond between them. We have found that this general model describes 
the order seen in LaMnOa for a range of parameters. In LaMnOs , as we explained above, 
the MnOg octahedra are not independent as they share oxygen atoms with their nearest 
neighbours (see Fig. El^a)). This introduces a constraint in the degrees of freedom for each 
octahedroni^. If a Mni-0 bond expands, the oxygen sharing 0-Mn2 bond should contract 
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and vice versa. According to that, Mn orbital ordering is not isotropic. Namely, if Mn-O-Mn 
was located along x direction in x — y plane, the 3(i orbit of Mn in the expanded Mni-0 
bond will be d^x3_^2 , x-state, but in the contracted 0-Mn2 bond it could be either d^y3_^2 
, y-state, or d^z^_,^2 , z-state, and the same order may exist in the next sublattice, as shown 
in Fig. n^b). The energy of a pair of identical orbits is Ji only if the lattice vector joining 
them is along the lobe of the orbit as shown in Fig. |21^b). The interaction energy of two 
identical orbits separated by a lattice vector that is not along the lobe is J2 as shown in 
Fig. m^c). In 3-state 3-dimensions anisotropic Potts model each site is occupied by one of 
the orbits x, y or z and the energy is described by the following Hamiltonian, 

1 ^ 

-y^aniso = - 2 5Z (/^ip'^^'^S, , (4) 

(id) 

where p = R- — R- is the direction of the nearest neighbour bond between two states. 

There are two exchange interactions for the anisotropic model. The 'head to head' 
interaction Ji is defined by Jx{p) = Ji for p = ±xa, Jy{p) = Ji for p = ±ya and Jz{p) = Ji 
for p = ±za. The 'side to side' interaction J2 is defined by Jx{p) = J2 ^01 p = ±ya or 
p = ±za, Jy{p) = J2 for p = ±xa or p = ±za and Jz{p) = -h for p = ±xa or p = ±ya. 
Thus each site has a coupling J2 to four neighbours and a coupling Ji to two neighbours. 
This is shown in Fig. Efb) and|2fc. It is worth mentioning that these types of interaction do 
not affect the overall cubic symmetry of the lattice. 

B. Monte Carlo Simulations 

This model has been studied using Monte Carlo (MC) simulations on 3-d finite lattices 
(with linear sizes L = 8, 10 and 12) for a range of values of Ji and J2 with periodic boundary 
conditions. All our simulations have made use of the Metropolis algorithm and with aver- 
aging performed over from 10^ to 10'^ Monte Carlo steps per site. Results were obtained by 
either cooling down from a high-temperature random configuration as discussed by Banavar 
et alM or heating up from the ground state. The results from the two procedures agree. 

Combination of the analytic methocP^ in the range Th < T < 00 and computational 
results in the range T < Th can be used to obtain the high-T limit to get the starting 
point, Tff, for our simulations. The high-T expansion of the internal energy as a function of 
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FIG. 3: Energy, u, per site versus (3\Ji\ (where (3 = l/ksT) for the anisotropic 3-state 3-d Potts 
model with L = 8 in the low P limit. The suitable starting point for our simulation is obtained, 
Th — 32\Ji\/kB- For simplicity we show the case where J2 = 0. 

temperature, u{T), from the Hamiltonian in Equation |3] for the case \Ji\ >> J2 is, 

u{T) = m(oo) - — ^ ( ^] , (5) 

^ ^ ^ ' 2kBTH V 4 y ' ^ ^ 

where z' is the number of the similar nearest neighbours, z' = 2, and u{oo) = '^^"^'^^ . The 
value of u{T) can be written as follow. 



«(T) 



Ji - 2 J2 A 
3 7^' 



(6) 



where A is a constant. Fig. IHl shows a plot of u{T) in units of Ji against For small 

P\Jil P\Ji\ = ii^< 0.03125, then, Th - 32\J,\/kB. 

We use the simulations to evaluate the thermodynamic quantities in the range < 
T < Th, and find the contribution for T > Th analytically. A full phase diagram with six 
regional phases in the Ji -J2 plane has been obtained numerically and analyticallyiSi, see 
Fig. H 



III. SIMULATION FOR LaMnOs 

The phase of interest in Fig. |3] is phase-5 where Ji < so that the head-to-head 
configuration is disfavoured and J2 > 0, so that, side-to-side configuration is favoured and 
< < 0.5. This has the correct low temperature pattern of orbital ordering seen in 
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FIG. 4: Ji - J2 phase diagraro=a of the orbital structures in a simple cubic lattice for general 
Hamiltonian (Eq. 0]). The ground state of LaMnOs occurs when Ji is AF, J2 is FM and < < 
0.5 which is phase-5 above. 

LaMnOs , as shown in Fig. ^b). In this phase we have an x,y checkerboard pattern in 
the X — y plane and this repeats itself so that the x and y states are in ordered chains up 
the z axis. The contribution to the energy comes from these ferromagnetic chains so the 
ground state energy is given by n(0) = — ^2- We find that at high-T the internal energy 
is u{oo) = |(Ji — 2J2) as expected from a random orbital array in the Potts model. The 
stabilisation energy of this phase is obtained as, 



The behaviour of the specific heat, Cy(T/|Ji|), in LaMnOs phase produced by this 
model (see Fig. Efa)) shows that it has a sharp peak corresponding to a first order transition 
in the range < J2/|-^i| < 0.5 and in the range < J2/\Ji\ < 0.25 there is a second, broader, 
peak. In the range 0.25 < < 0.5 the two peaks merge. The first peak obtained along 
the whole region of LaMnOa phase is at the transition from the orbital order to orbital 
disorder configuration occurring above Tjt- The second peak obtained along the region 
< J2/\Ji\ < 0.25 is caused by the short range order occurring because of the large value 
of the head-to-head interaction, Ji. We identify the short range order below the broad 



Au = u{oo) - n(0) = -(Ji - J2). 



(7) 
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FIG. 5: (a) Temperature dependence of specific heat, Cy, per site for the 3-state 3-d anisotropic 
Potts model at L = 12 where J2/\Ji \ = 0.01, 0.05, and 0.1 and Ji < 0. Tm is the temperature at the 
minimum value of Cy between the two peaks, (b) the dependence of the transition temperature 
on J2/I Ji|. 

peak with the short range order observed experimentally just above Tj^. In section ITVI we 
compare this short range order to the short range order at T = in the phase where = 0. 

We obtain the magnitudes of Ji and J2 as follows. The experimental value of Tjt is 
known, Tjt — 750K, and from Fig. Tjt = 750 = aJ2 where a is a constant which is 
obtained from the slope, J2 — 0.l2eV. We find Ji from the condition 0.05 < < 0.25, 
which gives Ji, lAeV > \Ji\ > O.SeV, which includes the value obtained by Millis et alM, 



The specific heat curves have a minimum at a temperature, T^, between the two peaks. 
We evaluate the entropy at by integrating the specific heat obtained from the simulation 
above and a temperature Th » Ji/ks where the entropy is given by kslogeS. 



The values of s(Tm) obtained in this way are plotted in Fig. IHlas a function of J2 / Ji- 

Figure El shows that ■s(Tm) values are almost independent of J2/\Ji\ along the two- 
peak region. The average value s{Tm) = (0.50 ± 0.02)kB is in good agreement with the 
experimental value obtained by Sanchez et al. who studied the local structure of LaMnOa 
across the JT transition at Tjt = 750K by means of X-ray absorption near edge structure 
and extended X-ray absorption fine structure at the Mn K-edge. They obtained from the 
heat capacity measurements of LaMnOa that the change in entropy between low T and 



O.QeV. 




(8) 
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FIG. 6: The entropy for LaMnOs phase versus J2/\Ji\ above Tjx where J2 = 0.0, 0.01, 0.05 and 
0.1 and Ji < 0. 

above Tj^, which is given as s{Tm) in our theory, is given by (0.515 ± 0.02)^;^. 

Since the relevant Potts model has J2/\Ji\ small it is instructive to consider, in next 
section, the limiting case of J-? = 0. 




0.2 0.4 0.6 0.8 
T/T 

JT 




0.2 0.4 0.6 0.8 1 
T/T 

JT 



FIG. 7: (a) The square of the order parameter, of the orbital ordering in LaMnOs vs T/Tjt 
at J2/\Ji\ = 0.05,0.10,15 and 0.20 with L = 12, (b) The experimental results for the square of 
the order parameter of the orbital ordering as function of the reduced temperature, T/Tjt, (after 
Murakami et al.^). 

The square of the order parameter of the orbital ordering, in LaMnOs is found 
as a function of the reduced temperature, from the simulation as shown in Fig. [7| 
We obtain the same value of ^"^{T/Tjx) from Monte Carlo simulations independent of the 
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direction of temperature change taking account of the fact that there are three nontrivial 
inequivalent ordered phases, corresponding to a choice of Q. 

For ordering characterised by Q = (110) we define a site (n,m,l) as belonging to the 
odd (even) sublattice if + m is an odd (even) integer. We define the probabihty, pi"'' , by 
defining the total number of odd sites that are occupied by an x orbital, A''^. 

= ^ (9) 

This definition is extended to define Py and p° and similarly for the even sites Py and pi- 
The order parameter of the orbital ordering for an x or y orbit in an odd site is. 



i"^ = p^:^-lip^;^+p^"^), (10) 



and 

1 
2 

Similarly, we are able to obtain the ordering on the even sites (when m + n is even) e.g. 



i°^=#-o(pi°^+Pi°^)- (11) 



<p^:^=pi'^-^ipi^+pi'^), (12) 

where p^x\py''' and pi^'* are the probability to find the even site occupied by x,y or z orbit 
respectively. The total order parameter of the orbital ordering in LaMnOa is, 

$(no)(T/r,^) = i^^L±M, (13) 

where (px = — 4>^x^ and (py = (p^y^ — (p^y \ We do not have a term for cpz because it is not 
included in $(iio)(T/T7j') when Q = (110). When we start the simulation from high-T, we 
find that the system goes into one of the three equivalent Q ordered phases, (110), (101) or 
(Oil) at the same value of Tjt- 

Figure [7(a) shows that the square of the order parameter of orbital ordering in this 
model phase has an unusual behaviour. The orbital ordering is essentially perfect for T < 
QATjT- It decreases through the range 0.4 < 7^ < 1 faster than an Ising model. Finally, 
it goes rapidly to zero at Tjt because there is phase transition at this temperature which 
appears to be first order (but this is not proved rigorously because of the finite size of 
the simulation). The order parameter is plotted as a function of the reduced temperature, 
T/TjT, in Fig. [7|^a). It may be seen that the order parameter is remarkably insensitive to the 
value of J2 / Ji provided that J2/J1 << 1- Figure [7||^b) shows the experimental results from 
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Murakami et al.^. There is a qualitative agreement between the simple model proposed here 
and the experiments. The size of the first order jump at Tjt is well reproduced by the Potts 
model. The experimental results have a strong concave region for 0.25 < 5^ < 0.9 which 
is reproduced qualitatively by the theory. This is unusual behaviour of an order parameter 
as in most cases (for example Ising model) the value of |^| increases monotonically with 
temperature. 

IV. PHASE FOR J2 = 

In this section we show how a study of the phase obtained for J2/\Ji \ = gives further 
understanding of the orbital ordering in LaMnOs just above Tjt- When we set Ji as AF 
and put J2 = (at J2 = there is no phase transition) a phase which has some physics of 
LaMnOs phase above Tjt is obtained. Namely, there is no longer side-to side interaction 
but the head-to-head configuration is still disfavoured. 

The energy varies from the high-T limit, u{oo) = Ji/3 to zero at T = where there 
are no head-to-head configurations. In comparison of the new model with AF Potts model, 
because the anisotropic case has fewer antiferromagnetic interaction we expect higher value 
for the ground state entropy, s(0), at T = 0. 

Therefore, the configuration of LaMnOs above Tjt is affected by the short range order 
obtained by the head-to-head interaction. It is known that the entropy of the 3-state Potts 
model is kslogS at high-T and nearly ^loge2 at T = 0.— The anisotropic Potts model has 
an entropy at T = equal to s(0) = (0.6 ± 0.02)/^^ which is dramatically higher than s{Tm) 
that for LaMnOs phase above Tjt through the range < ^ < 0.25. 

V. CONCLUSION 

A modified model called the anisotropic Potts model has been presented to study 
the orbital ordering and the local JT distortion above JT temperature, Tjt, in LaMnOs . 
The interactions between the three possible orbits depended both on the type of orbits and 
their relative direction on the lattice. Suitable values of the exchange interaction have been 
chosen to give the correct low temperature ground state for the orbital ordering phase in 
LaMnOs. The phase obtained from these interaction values has a first order transition at 
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the JT temperature, Tjt- 

Short range order was obtained above Tjt for a wide ratio of the exchange interac- 
tions < J2/\Ji\ < 0.25. The short range order reduced the entropy from the value of 
three randomly occupied orbitals, ksloQe'i = l.OdSks, to be (0.5 ± 0.02)kB which is in fair 
agreement with the experimental value obtained by Sanchez et alM. When the side-to-side 
interaction was put to zero, J2 = 0, we obtained a phase that does not order down to the 
lowest temperature. This phase had an entropy at low temperature that was higher than 
that obtained for LaMnOs phase above Tj^. 

We obtained the order parameter of the orbital ordering phases as a function of tem- 
perature in good agreement with the published results. The square of the order parameter, 
^"^{T/Tjt), of the orbital ordering decreases faster for T > OA /Tjt than an Ising model. 
At Tjt it goes to zero discontinuously in a first order transition. 

The ratio of J2 /Ji necessary to give the LaMnOs ground state are < J2/\Ji\ < 0.25 
and within this range we find that the results are almost independent of the exact value of 
J2 /\Ji\- In particular we find similar shapes of the curve $^ as a function of the reduced 
temperature, T/Tjt, as shown in Fig. Eland also that entropy above Tjt as shown in Fig. IHl 
is also approximately constant. Furthermore, the range of values of J2 /\Ji \ required to give 
the observed ordering temperature is consistent with those proposed from a microscopic 
model. Thus we have shown that the temperature dependence of the orbital ordering in 
LaMnOs is represented well by a phenomenological three state Potts model. 
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